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An approach to experimentally exploring electronic correlation functions in mesoscopic regimes 
is proposed. The idea is to monitor the mesoscopic fluctuations of a tunneling current flowing 
between the two layers of a semiconductor double-quantum-well structure. From the dependence 
of these fluctuations on external parameters, such as in-plane or perpendicular magnetic fields, 
external bias voltages, etc., the temporal and spatial dependence of various prominent correlation 
functions of mesoscopic physics can be determined. Due to the absence of spatially localized external 
probes, the method provides a way to explore the interplay of interaction and localization effects in 
two-dimensional systems within a relatively unperturbed environment. We describe the theoretical 
background of the approach and quantitatively discuss the behavior of the current fluctuations in 
diffusive and ergodic regimes. The influence of both various interaction mechanisms and localization 
effects on the current is discussed. Finally a proposal is made on how, at least in principle, the 
method may be used to experimentally determine the relevant critical exponents of localization- 
delocalization transitions. 

PACS numbers; 73.21.-b, 73.23.Hk, 73.40. Gk, 73.50.-h 



I. INTRODUCTION 

At low temperatures, disordered or chaotic electronic 
systems are strongly affected by mechanisms of quantum 
interference. These interference effects manifest them- 
selves in anomalously strong fluctuations of both ther- 
modynamic and transport observables and in the spatial 
localization of quantum mechanical wave functions [Q. 
They find their common origin in an interplay of the 
classical nonintegrability of the charge carrier dynamics 
and the wave nature of quantum mechanical propagation. 
To quantitatively characterize this 'mesoscopic' behavior, 
one commonly employs correlation functions of the type 

i^(xi,X2,X3,X4;w, Aa) = 

= (G" (xi, X2; e, a)G+(x3, X4; e -I- w, a + Aa)) , 

where is the retarded/advanced single-particle Green 
function. Here (. . .) stands for some kind of averag- 
ing (e.g., averaging over realizations of disorder) and the 
parameter a symbolically represents an optional depen- 
dence of the Green function on external control param- 
eters (magnetic fields, gate voltages or others). Corre- 
lation functions of this type appear as the 'most micro- 
scopic' building block in the analysis of the majority of 
fluctuating mesoscopic observables. As a consequence of 
constructive quantum interference these objects become 
long ranged whenever the spatial arguments Xi , . . . , X4 
are pairwise close (on scales of /min, the range of the av- 
eraged Green functions {G"^)). Specifically, 

• For X = xi w X2 and x' = X3 w X4, F['^1(x, x') = 
F(x, X, x',x') describes the fluctuations of the (lo- 
cal) density of states (DoS), and, thus, the thermo- 
dynamic fluctuations and parametric correlations. 



• For Aa = 0, i^[^l(x,x') = F(x,x',x',x) describes 
the total probability of propagation from x to x'. 
This is the generalized 'diffuson', a quantity of key 
relevance in the context of mesoscopic transport. 

• Finally, in a system with unbroken time rever- 
sal symmetry, F['-^](x, x') = i^(x, x', x, x'), the 
cooperon, becomes long ranged, too. 

The dependence of the correlators F^'*'^''-^] on the long- 
ranged distance r = |x — x'|, the energy difference uj, 
and Aa fundamentally characterizes a multitude of meso- 
scopic phenomena For this reason, many theoretical 
investigations in mesoscopic physics concentrate on an 
analysis of these objects. Experimentally, however, it 
has proven difhcult to access the correlation functions 
p[d;D;C] directly: Ideally, one would like to continuously 
measure the dependence of the correlators F^''''^'^^ over a 
range of at least the parameters r and u. Irritatingly, this 
cannot be achieved within experimental setups based on 
a standard device-contact-electron system architecture. 
In fact, the mere presence of local contacts introduces 
an entire spectrum of difficulties obstructing the contin- 
uous experimental spectroscopy of transport and spectral 
correlation functions: First, the fixed attachment of lo- 
cal current/ voltage electrodes prevents one from contin- 
uously monitoring the scale (r) dependence of transport 
correlation functions. This problem does not exist in 
measurements based on local tunneling tips In those, 
however, the electronic state of the tunneling device as 
well as its coupling to the electron system have to be 
precisely known to draw quantitatively reliable conclu- 
sions on the nature of the bulk electronic correlations of 
the latter. In particular, the distance between the device 
and the electron system has to be kept constant with 
atomic precision. These conditions can hardly be met 
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under realistic conditions. Second, both local contacts 
and tunneling tips tend to disturb the electron system 
under investigation. In interacting systems, they lead to 
various manifestations of the orthogonality catastrophe. 
As a consequence, much of the measured current/voltage 
characteristics describes the process of local accommo- 
dation of charge carriers at the interface, rather than 
the electronic correlations of the bulk system. Third, 
a division between system and contacts of a mesoscopic 
conductor is, to a large extent, arbitrary. Quantum in- 
terference phenomena in mesoscopic systems tend to be 
highly nonlocal in space, and often it is not clear, where 
the physical processes responsible for the outcome of an 
experiment took place, in the 'device', the 'contacts', or 
all over the place. 




FIG. 1. Schematic setup of the approach: two parallel 
two-dimensional electron systems are subject to a bias volt- 
age and external magnetic fields. A tunneling current flowing 
between the layers is recorded as a function of the external 
control parameters. 

These problems with contacted systems and local tun- 
neling devices led to the idea to use a second electron 
system with essentially known properties as an extended 
tunneling spectroscope Setups of such type can 

conveniently be realized in double-quantum-well struc- 
tures embedded in semiconductor heterostructures Q , as 
shown schematically in Fig. [l]: Two electron systems con- 
fined in parallel wells separated by an isolating barrier of 
uniform thickness d ^ 10—20 nm form a double-layer 
system of two two-dimensional electron systems. A tun- 
neling current / from one layer to the other is driven by 
applying a voltage difference V between them. The tun- 
neling region is typically a few 100 x 100 iiir? in extent. 
In the absence of any disorder scattering and/or tunnel- 
ing amplitude inhomogeneities, the tunneling from one 
layer to the other can only occur if energy and momen- 
tum are conserved. This leads to the resonant behav- 
ior of the tunneling current that is characteristic of two- 
dimensional systems. More generally, for constrained ge- 
ometries (e.g., quantum wire/two-dimensional electron 
system) the field dependence of these resonances can be 
analyzed to obtain information on the dispersion of the 
fundamental excitations in the two systems |^ . 

However, in 'real' systems, inhomogeneities in the tun- 
neling barrier thickness, static disorder, and other non- 
momentum conserving imperfections will lead to mod- 
ifications of the idealized resonant current profile. Of 



these intruding mechanisms, the first appears to be the 
most serious: the current will respond with exponential 
sensitivity to any fluctuations of the layer separation; 
for strong enough spatial variations one may run into 
a scenario, where tunneling occurs only at a sparse set 
of 'hot spots', with no traces of a resonant profile left 
. (Some characteristics of this type of current fiow will 
be discussed below.) However, recent technological ad- 
vances have made it possible to manufacture double-well 
systems with near-monolayer precision. In such devices, 
fluctuations in the tunneling matrix elements are reduced 
down to values of 0(10%) and can be absorbed into 
a renormalization of the effective in-plane disorder. In 
the present paper, the focus will be on transport in these 
near-planar devices. 

Even if the tunneling is homogeneous, static disorder 
will broaden the resonant behavior and introduce fluctu- 
ations. The broadening of the average current is related 
to the dynamics on short time scales In contrast, 
the fluctuations contain information about physical pro- 
cesses on much larger time scales of the order of, e.g., 
the diffusion time through the system. It is the purpose 
of this paper to investigate the nature of these fluctua- 
tions and their relation to the aforementioned electronic 
correlation functions. 

In fact, we will see below that detailed information on 
the correlation functions F^'*'^''^! can be extracted from 
the tunneling current fluctuations without disturbing the 
system. Moreover, (i) the tunneling takes place uniformly 
at all points of the layers which means that an averaging 
over spatially fluctuating structures (e.g., details of the 
microscopic wave function amplitudes) is intrinsic to the 
data contained in the current, (ii) Several parameters can 
be tuned to gain information: The bias voltage resolves 
energetic correlations, a parallel magnetic field resolves 
spatial correlations, and a perpendicular magnetic field 
may serve as a control parameter for parametric corre- 
lations (i.e., correlations between Green functions evalu- 
ated at different values of external control parameters), 
(iii) The geometry of the layers can be designed freely, so 
that it is possible to study different regimes of particle 
dynamics (e.g., ergodic, ballistic, diffusive, etc.). 

In this work, after the introduction of the general theo- 
retical background of the tunneling current statistics, we 
will consider two prototypical system classes. First, we 
study extended systems for which the phase coherence 
length is much larger than the microscopic length 
Irain and a crossover or a transition from diffusive motion 
to Anderson localization may take place. For such sys- 
tems, a parallel magnetic field By can be employed as an 
instrument for resolving the long-range behavior of the 
correlation functions fI^'*-^!. Importantly, the field align- 
ment parallel to the two-dimensional planes implies that 
the charge carrier dynamics is not affected by B||. While 
in the diffusive regime explicit expressions for the correla- 
tion functions are known, no quantitative expressions for 
regimes with strong (nonperturbative) localization and 
interaction effects are available. However, in a regime of 
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localization-delocalization crossover or transition, scal- 
ing behavior is expected to restrict the functional form 
of fI^'^^I. This opens the possibility to extract the rel- 
evant critical indices from tunneling conductance mea- 
surements. 

Second, we study a geometry, where one of the layers 
forms a ballistic quantum dot (L < /min) in the ergodic 
regime. The other, extended, layer serves as the spec- 
trometer. For the quantum dot, the parametric corre- 
lations with respect to a perpendicular magnetic field, 
present in F^'^'^''-^\ can be obtained from the current 
fluctuations. A similar setup has already been realized 
experimentally by Sivan et al. |p!o| . In that work, a sin- 
gle level (in contrast to our extended system) was used 
as a spectrometer to study a quantum dot device. This 
experiment led to results for the functional form of the 
correlator F^'^\ compatible with theoretical predictions 
from random matrix theory. However, one would expect 
that the data obtained from single-level spectroscopy 
is still weighted with nonuniversal wave-function ampli- 
tudes specific to the isolated 'monitor level'. In contrast, 
for the two-dimensional layer/quantum dot setup con- 
sidered here, the current flow is extended and spatially 
uniform. As a consequence, the tunneling current fluctu- 
ations are microscopically related to the purely spectral 
content of parametric correlations. Below we will estab- 
lish the quantitative connection between the field and 
voltage dependence of the tunneling current fluctuations 
and a number of correlation functions that have been an- 
alyzed in the recent theoretical literature ||ll|. Moreover, 
we will try to assess to what extent these connections, ob- 
tained for the chaotic noninteracting electron gas, may be 
susceptible to interaction mechanisms such as Coulomb 
drag 1 12 1^ or Coulomb blockade effects Q. 

The outline of the paper is as follows: In Sec. || we 
specify our model and review the general formula for the 
tunneling current. Turning to new results, we show that 
the Fourier transform of the tunneling conductance corre- 
lator with respect to the magnetic field directly yields the 
spatially resolved correlation functions F^^''-^^ In Sec. |l| 
the theory will be applied to diffusive and anomalous dif- 
fusive systems, and in Sec. IV to an analysis of spectral 
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and parametric correlations in finite quantum dots. The 
impact of Coulomb charging effects on the tunneling cur- 
rent fluctuations will be discussed in Sec. 
elude in Sec. 



II. THEORY OF TUNNELING CURRENTS 



A. The current formula 



Consider a double-layer system consisting of two par- 
allel two-dimensional electron gases (2DEGs), labelled 1 
and 2, respectively (see Fig. ||). The two layers are sep- 
arated by a tunneling barrier that we assume to be uni- 
form. We aim to analyze the tunneling current / under 



conditions, where the tunneling is weak (in a sense to 
be specified momentarily). After matching the electron 
densities in both layers by adjusting their chemical po- 
tential /i, the current becomes a function of bias voltage 
V, temperature T, and, optionally, a magnetic field B, 
I^I{V,B;T,^i). 

Quantum mechanically, the system can be represented 
in terms of a tunneling Hamiltonian, H = Hi + H2 + Hq-, 
where Hi and H2 describe layer 1 and 2, respectively, 
while Hj- describes the transfer of electrons between the 
layers . Choosing a gauge where the bias voltage has 
been transferred to the tunneling matrix elements, Hq- 
can be written as 



IT 



-iVt 



V;\(x)^2(x')+h.c.) , (1) 



where 7^x' is the tunneling amplitude from x in layer 1 to 
x' in layer 2, and ipj are electron creation and anni- 
hilation fields for layer j = 1, 2. For convenience, we use 
units where h^c = kB^e — l- Since the tunneling matrix 
elements 7^x' decrease exponentially (on atomic scales) 
as a function of |x — x'| we model T^x' as a spatially local 
object, T^x' = ^<5(x — x'). In the weak tunneling regime, 
i.e., when the typical time after which an electron tun- 
nels is larger than the characteristic time scale that is 
to be resolved in the experiment, single-tunneling events 
dominate |l6| . Then, the tunneling current reads ||l^,^ 

liV, B) = 2 y" d^x (fx' J (del^l) TxTx*, e^^^''-''') x 

X Ai(x, x'; e, B^)A2{x', x; e - B^) , (2) 

with the abbreviation /(del^l) = /^||K(e ~ V) ~ 

npie)]. Here npie) = (l + e^"""^)/^)"^ is the Fermi dis- 
tribution function at temperature T and Fermi energy 
ep. The characteristic momentum scale set by the par- 
allel field is q = c?B|| x e^, where is the unit vector 
perpendicular to the plane. 

In Eq. (^) the quantities of main interest are the spec- 
tral functions of layer j, Aj = i{G'^ — Gj), where 
is the retarded/ advanced one-particle Green function. 
These spectral functions depend on the bias voltage V 
applied to layer 2, and on the perpendicular magnetic 
field B±. In fact, it will be our main objective to ob- 
tain information on these objects through their parame- 
ter dependences. In this context, it is crucial to note that 
B|| does not change the dynamics within the individual 
layers, but merely weighs the tunneling current with an 
Aharonov-Bohm-type phase. The sensitivity of the cur- 
rent to this flux will help to gain information about the 
long-range propagation within the layers. A caricature of 
the basic idea is depicted in Fig. |l|. This figure illustrates 
the basic physical processes underlying the current flow 
as described by Eq. (||): An electron tunnels at point 
x from layer 2 to 1, leaving a hole behind. It propa- 
gates within that layer to point x', where it tunnels back 
to layer 2 and recombines with the hole. The in-plane 
magnetic field enters the formula via the flux through 
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this electron-hole loop. Therefore, the in-plane magnetic 
field dependence of the tunneling current contains infor- 
mation about the typical area enclosed in the loop that 
in turn is determined by the typical range of propagation 
within the layers. 

To further simplify the analysis, we note that in the 
absence of significant interaction corrections the spectral 
functions themselves do not exhibit temperature depen- 
dence. Under these conditions, a simple integral relation 
between currents at zero and finite temperatures holds: 



^)/(T = 0,eF = e). (3) 



/(r,eF) = - de 



In the following, unless stated otherwise, all results will 
be given for T — only. The generalization to finite tem- 
perature - essentially a smearing of the T — results - 
obtains from Eq. (||). 

In this paper, we are primarily interested in the tun- 
neling current flowing between disordered systems. The 
microscopic properties of these systems will be described 
by some disorder distribution function about which we 
make three idealizing assumptions: 

(1) The disorder potentials of the two layers are essen- 
tially uncorrelated fig] . 

(2) The e-e interaction and higher order tunneling pro- 
cesses are not able to introduce significant interlayer cor- 
relations in the motion of the charge carriers. Practi- 
cally, this means that impurity averages can be taken for 
each layer independently. Roughly speaking, e-e interac- 
tion effects can be divided into three groups: momentum 
transfer between the layers ('Coulomb drag'), charging 
effects associated with the tunneling process, and self- 
energy corrections (due to inelastic scattering and de- 
phasing). Sizable interlayer Coulomb correlations may 
arise in very clean systems subject to strong perpendic- 
ular magnetic fields. In such systems, the e-e interac- 
tion can stabilize a fractional quantum Hall phase p9| ] 
and interlayer e-e interactions lead to additional corre- 
lation phenomena (spontaneous coherence and quantum 
Hall ferromagnets ||2^). In contrast, for stron g en ough 
disorder, the e-e interaction is less significant and, 
thus, the effective random potentials in each layer can be 
treated as statistically independent of each other. This 
assumption can be tested experimentally, as will be dis- 
cussed below. Charging effects will inevitably influence 
the tunneling current at low bias; we p ostpone a discus- 
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sion of these corrections until Sec. 

(3) Disorder does not significantly affect the spatial 
homogeneity of the tunneling, i.e. the tunneling ma- 
trix elements do not depend on position, 7^ = T. The 
spatial homogeneity of the tunneling is very sensitive to 
the thickness of the barrier. However, it is now possi- 
ble to grow heterostructures with near-monolayer preci- 
sion and, thus, achieve a tunneling probability that is 
almost spatially constant p2p^ ]. (In high precision de- 
vices, the space dependent relative fluctuations in the 
tunneling probability can be reduced to about 10% and 



lower 0.) The validity of this assumption can be tested 
experimentally. In the extreme case, where tunneling oc- 
curs only through tunneling centers or 'pinholes' Q, the 
resonant behavior of the tunneling current disappears. If 
the spacing between pinholes exceeds the mean free path, 
the average current is just proportional to the product of 
the local densities of states, (/) ~ vii'2, and, furthermore, 
becomes independent of magnetic field. 

In principle, a full (angle-resolved) analysis of the 
field dependence of the current would allow one to ex- 
tract information about the distribution of tunneling cen- 
ters [|4j. However, further discussion of this type of 'tun- 
neling center spectroscopy' is beyond the scope of this 
paper. Keeping in mind that the working assumption of 
near-homogeneous tunneling can be put to experimen- 
tal test, we hereafter concentrate on the case 7^ = T. 
(Note that, although significant inhomogeneities would 
largely obstruct the detection of transport correlations, 
they only have a minor effect on the analysis of spectral 
correlations.) 



B. Average current 

To foster the discussion of fluctuations below, this sec- 
tion briefly recapitulates the main characteristics of the 
average tunnel current p5[ . Under the assumptions for- 
mulated above it is given by 



{liV, B)) = ^ \r\^ jde j d\ e^i" 



X (Ai(x; e, Si))(A2(-x; e - V, B^)) . (4) 

The current, Eq. (^, is characterized by the averaged 
one-particle Green function (G'*j(x,x'; e)). This quan- 
tity is short-ranged on a For small perpen- 
dicular magnetic fields B±, Zmin is set by the mean free 
path / = vpT^ where vf is the Fermi velocity and t the 
elastic scattering time. (In cases, where the scattering 
times in the layers are different, we need to generalize to 
Tj, j = 1,2.) We here focus on systems, where the mean 
free path I is much shorter than the linear system size 
L |26[ |. For stronger perpendicular magnetic fields, with 
cyclotron frequency lOc = Bj_/m exceeding the inverse 
elastic scattering rate r~^, the classical cyclotron radius 
Rc = vf / is smaller than the mean free path / implying 
that ?min = Rc (see, e.g., Ref. 

The average current has been studied theoretically ||] 
and experimentally |2^,^. For matched Fermi energies 
and vanishing magnetic fields, the theoretical result reads 



a(^,B^O)) 
V 



Go 



r2 



Go 



(5) 



where F = Fi -h Fj, F^- = 1/(2tj) is the hne-width of the 
Lorentzian-shaped average spectral function in layer j, v 
is the single-particle-level density of states per unit area. 
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and Go is the characteristic low-bias average conductance 
of the system. For the average current, the condition of 
weak tunneling is 



Tr <r, 



(6) 



where Tr = \T\ is the inverse of the (golden rule) 
rate at which a particle propagating in one layer tunnels 
to the other. Henceforth we will focus on the regime 
of small bias voltage V ^ T. Under these conditions, 
(/(y, B = 0))/V « Go. In the presence of a moder- 
ately weak in-plane magnetic field (|q|~"'^ ^ Fermi wave- 
length), this generalizes to 



a(^,B)) 

V 



= Go/(|q|0, 



(7) 



where the scaling function / exhibits the asymptotic be- 



havior f{x < 1) = l+0{x^) and f{x > 1) 



In the 



following, we will concentrate on the weak field regime, 
|q|Z«l^/«l. 

Before moving on to our main issue, mesoscopic fluctu- 
ations of I{V, B), let us make a few more remarks on the 
average current. Specifically, we wish to argue that from 
the known behavior of (/) , conclusions on the validity of 
some of the assumptions made above can be drawn. In 
Refs. [|2| and |Q the influence of a perpendicular mag- 
netic field on the average current between high mobility 
samples was investigated. It was found that strong fields 
lead to a suppression of the differential tunneling con- 
ductance, Gt = dl/dV, at zero bias. This phenomenon 
is called the 'tunncling-gap'. The splitting of the con- 
ductance peak a,t V — 0, characterized by some field- 
dependent gap energy A{B), is due to the Coulomb in- 
teraction within and between the layers. As shown theo- 
retically in and confirmed experimentally in Refs. ||2^ ] 
and , the total current flowing at the split peaks (po- 
sitioned at = i^) equals the peak current at = in 
the absence of interactions. This observation implies that 
interactions in these experiments largely manifest them- 
selves in the form of self-energy corrections to the single- 
particle poles. In contrast, if strong interlayer correlation 
effects were present, the peak current would increase. In- 
deed, for strong interlayer correlations, the momenta of 
the particle and the hole constituting the 'current-loop' 
would be partially correlated. This should lead to grad- 
ual resurrection of the resonant behavior characteristic 
for the clean case and, therefore, to an unsplit zero-bias 
conductance peak. 

Even at B = the diffusive zero-bias anomaly leads to 
some splitting of the peak. However, the separation of the 
maxima of the subpeaks is small in 1/ y^g ||2^ and will be 
neglected here. Note that, when lowering ff, nonpertur- 
bative (in 1 /g) zero-bias effects can occur pQ]. Finally, 
in Sec. [VB| the analog of the zero-bias anomaly in finite 
quantum dots (Coulomb blockade) will be considered. 



C. Fluctuations 

We next turn our attention to the fluctuations of the 
tunneling current. As a starting point we will use the 
formula 

I{V,B) = ^\T\^ Jd^xd^x' J dee*'i("-"') X 

X Ai(x, x'; e, B^)A2{^', x; e - B^) (8) 

for the zero-temperature current at uniform tunneling 
probability. As we are interested in correlations on large 
time scales, Eq. (||) for the range of applicability of the 
weak tunneling approximation has to be replaced by the 
more restrictive condition 



Ft < max(r, ^,V,VF\q\) 



(9) 



where is the phase coherence time and v-p\q\ the char- 
acteristic energy scale set by the parallel magnetic flcld. 
This inequality states that the probability for an electron 
to tunnel, while moving coherently within one layer, is 
low. 

To describe fluctuations of the current and related 
quantities, we will consider correlation functions of the 
type 



Cx{z,z') 



{X{z)X{z')), 



where X is an observable, {AB)c = {AB) — {A){B), and 

z''^ represents the set of parameters {V'^'\ B^'^ ep^}. [To 
avoid confusion, let us reiterate that in this paper an- 
gular brackets stand for averaging over an external set 
of parameters, not for a quantum mechanical average. 
For example, in our discussion below, X{z) = Gr(V^, B) 
may stand for the conductance measured at a certain 
field/voltage configuration. The subsequent (. . .)-average 
will then be over configurational fluctuations.] The sup- 
pression of the parameter dependence in the normaliza- 
tion denominator indicates that, on the ^-scales relevant 
for the structure of fluctuations, the parameter depen- 
dence of the averaged observables is negligibly weak. 

In most of the following, we will concentrate on the 
correlation function Co of the differential conductance 
Gt. This quantity is (a) experimentally more relevant 
than the current correlation function G/ and (b) tends 
to exhibit more pronounced structure. Indeed, it is 
straightforward to show that the current and the con- 
ductance correlation function, respectively, are related 
through Ci{V) = 2V-^ doj {V - w) Gg(w), i.e. G/ is 
obtained from Cg through an integral average. 

Averaging (the square of) Eq. (ra) over disorder, one 
verifies that 



Gg(c^, B,B') 



^{Gt{V--B)Gt{V'-,-B')) 



5 



— ^I^^l f TT ^2 iq(xi-X2)+iq'(x3-X4) 

x3fii^i(xi,X2,X3,X4;2)5ftF2(x2,xi,X4,X3;z) , (10) 

where Go = (Gt) as discussed in the previous section. 
The objects 

-F'l/2(Xl,X2,X3,X4; z) = 

where z — {w, Bj_, B'j^}, and uj — V — V, are the basic 
two-particle correlation functions discussed in the Intro- 
duction. That Eq. ( ]To| ) contains the product of two of 
these correlators is a direct consequence of our assump- 
tion of negligible interlayer disorder correlation. Notice 
that while the correlation functions F of noninteracting 
systems categorically depend only on the energy differ- 
ence between the two Green functions, the dependence 
on the perpendicular fields can be more complicated. 




FIG. 2. Diagrams contributing to Cg//: a) density-density, 
b) diffuson, and c) cooperon. The thick wavy lines denote the 
diffuson and cooperon propagators. 

The fourfold integration over the coordinates Xi 
implies that all three contributions discussed above, 
density-density F^fK diffuson F^^\ and cooperon F^'-^\ 
contribute to Eq. ( [l0| ) (see Fig. ^). At this stage, the role 
of the weak in-plane magnetic field becomes clear. As dis- 
cussed above, the correlators F^''-'^'^^ are long ranged (as 
compared to the microscopic spatial extent of the average 
Green functions contributing to (/)). This means that 
Eq. (|l^ is field sensitive - through the magnetic wave 
vector - on small magnetic field scales. The characteris- 
tic field strength is determined through |q| = dB^^ ~ Lj^, 
where Li^ is the typical distance a particle propagates 
during time For example, for a medium charac- 

terized by diffusive motion with diffusion constant D, 
L^j (DjiSf-^"^ . Using that for the three fundamental 
correlators the coordinates are pairwise equal [with an ac- 
curacy of ©(/mill)] and neglecting factors ^ |q|^min ^ 1, 
Eq. ( p^ ) assumes the form 



Cg(c^,B,B') = 




X (^S^-'^)^^-^'^ sRi^I^l (x, x'; z) 3fii^i^l (x', x; z) + 

^g»(q+q')(x-x') gfiF[t^l(x,x';z)3fii^f l(x',x;z) + 

-f5RF|''l(x,x';z)5RFi''l(x',x;z)) , (11) 

where g — I-kDv is the dimcnsionlcss conductance. 
Equation ( [ll| ) states that the diffuson contribution 
couples to the difference, B|7 = B|| ^BJ,, of the two 



in-plane field vectors, the cooperon contribution FI*^! to 
the sum, Bjj" = B|| + B'|, whereas the density-density 

contribution F^"^^ is B|| -insensitive. (Later on we will see 
how information on F^'^^ can be extracted from the de- 
pendence on the perpendicular field B\^.) Equation (|ll| ) 
holds true for extended systems, where the unconstrained 
integration over x, x' implies momentum conservation, as 
well as for restricted systems, where the in-plane momen- 
tum is not conserved in tunneling. 

Finally, if both systems are extended, Fourier trans- 
forming Eq. ( pd] ) in the magnetic field we obtain 

5i^^[-l(Ax;.)5R^^fl(-Ax;z)=(^)'x 
y d2B+e-*''(^"^^^^GG(c^,B+,B-), 

^Ff\^^-z)^Ff\-^^-z) = (^)' X 

y rf2B-e-''^(^-xBp^c'G(a;,B+,B-), (12) 

where Ax is the difference between the two spatial argu- 
ments of the correlation functions F, and we have used 
that qx = (i(x x Bii)^ as well as the result (^) for Gq. 

Equation ( p^ ) contains a central message of our paper: 
Detailed spectral and spatial information on the correla- 
tion functions F can be obtained from the dependence 
of the tunneling current on a parallel magnetic field. (In 
contrast to contact measurements,) the current approach 
to exploring correlation functions enables one to contin- 
uously measure spatial scale dependences, and does not 
incorporate strong local perturbations. If one of the lay- 
ers is a finite quantum dot, Eq. ( pi] ) gives the general 
relation between the current fluctuations and the spec- 
tral correlation functions. In the next two sections, we 
will discuss applications of this general concept to some 
concrete problems. 



III. ANOMALOUS DIFFUSION 

In this section, we are going to apply Eq. ( |l^ ) to the 
problem of (anomalous) diffusion in spatially extended 
structures. We first note that for the limiting cases of 
purely ballistic and diffusive dynamics, respectively, the 
correlation functions fI-^'*-^! can be calculated explicitly. 
For ballistic systems, a straightforward integration over 
the momenta of the single particle Green functions ob- 
tains 

5iF[^^^l(r;c^)~— cos— . (13) 

For diffusive systems, leading-order diagrammatic per- 
turbation theory (one diffuson/cooperon approximation) 
leads to 
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5Rf^[^'^l(r;tj)^ — kerf-^^ , (14) 

where ker(x) is the Thomson function ]3l]] . For small 
x, this function can be approximated as ker(x) ~ — C — 
ln(x/2) (C « 0.58 Euler's constant). The ellipses stand 
for weak-localization-type contributions of higher order 
in the number of diffusons and cooperons. These cor- 
rections scale with negative powers of the dimensionless 
conductance g. By definition, we will denote a system as 
'diffusive' if g 3> 1 and weak localization does not play a 
role. 

To get some idea about the strength of the tunneling 
current fluctuations let us briefly discuss the B = cur- 
rent correlation function Ci{V) for two different setups: 
(a) two disordered layers and (b) only one disordered 
layer and one 'clean' layer. Here 'clean' means that / is 
much larger than system size L. Substituting Eqs. ( p^ ) 
and (HJ) into Eq. (|lO|) and integrating over frequencies 
one finds in case (a) 

Ci{V) = g-^^ic^"^ + 2c[^l HV/En.)] , (15) 

where ETh = gS is the Thouless energy. Eq. ([l^) has 
been derived under the assumption V > E^h- Physi- 
cally, this means that on time scales t ~ V~^, the charge 
carriers do not have enough time to diffusively explore 
the entire system area. For smaller voltages, a crossover 
to an ergodic regime, discussed in the next section, takes 
place. The two numerical coefficients d'^^ = 9/(87r) and 
cl^l — l/(47r) determine the strength of the density- 
density and diffuson contribution, respectively. The fac- 
tor of 2 multiplying c'^l expresses the fact that in the 
field-free case, the diffuson and cooperon contribution, 
respectively, are equal and add. For case (b), the ex- 
pression lo oks s imilar, however, instead of the logarithm 
a factor \fV[V appears. This means that in the regime 
of interest, ^ <C F, the current fluctuations between a 
clean and a disordered layer are largely due to fluctua- 
tions in the density of states. This result can easily be 
understood qualitatively: In a clean system, the charge 
carriers move much faster than in a disordered system. 
As a result, the particles propagating in the disordered 
system do not have enough time to diffusively travel over 
large distances. This in turn implies that the diffuson and 
cooperon contribution to the correlation are reduced by a 
phase space reduction factor, whereas the density-density 
contribution, involving only Green functions taken at co- 
inciding points (within the clean system), remains un- 
affected. Actually, the density-density contribution to 
the current correlation is proportional to the variance of 
the number of levels in an energy window of V . This is 
very similar to the conductance fluctuations in conven- 
tional transport, which are related to the level number 
variance in an energy window of the size of ii'Th- For 
completeness, we mention that for case (b), c^'^l 4/7r 
and cl-^l = 16/(37r). Note that the conductance is self- 
averaging (~ i^^) in the thermodynamic limit. Further- 



more, the fluctuations are suppressed by a factor g~'^ . 
This is a phase-space reduction factor expressing the fact 
that to obtain averaging-insensitive contributions, two of 
the four spatial arguments of the correlation function F 
must be close to each other on scales I, cf. Eq. ([Tl|). 
Finally, we notice that already weak perpendicular mag- 
netic fields of 0(1/Aa;^), where Ax^ is the characteristic 
area of extent of fI'^I , suffice to suppress the cooperon 
contribution. This means that the dependence of the 
current fluctuations on a perpendicular field can be used 
to determine the maximum range of the correlation func- 
tions at frequencies u; ^ V. For V < t^^ , this scale is set 
by the dephasing length, L^. In analogy to the classical 
experiments by Bergmann [Q, the field dependence of 
the current for these low voltages can be used to estimate 

What can be said about systems with more com- 
plex types of dynamics, i.e., systems where localization 
and/or interaction corrections play some role. In prin- 
ciple, both weak localization and interaction corrections 
can be taken into account perturbatively, where the in- 
verse of the dimensionless conductance, g^^, represents 
the expansion parameter p9t . As g is lowered, these 
nondiffusive corrections become stronger and eventually, 
for g — 0(1), the perturbative description breaks down. 
However, relying on concepts of scaling theory, it is still 
possible to make some general statements about the be- 
havior of the strongly disordered electron gas: For 5^1, 
localization phenomena begin to qualitatively affect the 
dynamics. According to the one-parameter scaling the- 
ory |3^, the weakly interacting electron gas eventually 
flows into a localized regime provided that (a) spin-orbit 
scattering is negligible and (b) no strong perpendicular 
magnetic fields are present. In contrast, systems with 
significant spin-orbit scattering are expected to exhibit 
a true metal-insulator transition at some critical value 
gc ~ 0(1) [Q; 2D electrons (interacting as well as non- 
interacting) subject to strong magnetic fields undergo a 
metal-insulator transition responsible for the quantum 
Hall effect ]35[ |. Finally, in a number of experiments on 
2D electrons with strong interaction parameter > 1 
transport behavior has been observed that resembles a 
metal-insulator transition too. 

In all these phenomena (except, perhaps, the not suf- 
ficiently well-understood transport phenomena discussed 
in Ref. Pql) the concept of 'anomalous diffusion' plays a 
key role ||37|. Prior to the onset of strong localization, 
the electron dynamics undergoes a crossover from ordi- 
nary diffusive (5 2> 1) to anomalously diffusive {g ^ 1). 
Quite generally, the correlation function of anomalously 
diffusive electrons has the scaling form 

F[^l(r,a;)^(^£ye-2^/«, (16) 

where ^ is the localization length and rj a characteristic 
exponent related to the multifractal nature of states that 
are neither regularly extended nor fully localized 
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The length L^^ is related to the energy u by the so-called 
dynamical exponent z, 



-l/z 



(17) 



Whereas for noninteracting systems z = 2 as in ordinary 
diffusive systems, the value of z for interacting systems 
is controversial p9| . In systems with a true localization- 
delocalization transition, the localization length diverges 
with a characteristic exponent v upon approaching the 
transition point. In Eq. ([l^), we have assumed that ^ is 
smaller than L^. In the opposite case, would be the 
scale of exponential decay of the correlation function. 

According to Eq. (p^), the 'non-diffusivity' of the elec- 
tron dynamics can be characterized in terms of the three 
exponents z, v, and rj. To obtain these quantities one 
needs to know both the spatial and the energetic profile 
of the correlation function. In fact, the aforementioned 
difficulty to continuously monitor the spatial structure 
of electronic correlations has prevented previous exper- 
iments from determining the exponent 77. In contrast, 
the basic relation ([l^ ) does, at least in principle, con- 
tain all the information needed to extract all exponents 
of anomalous diffusion. In the following, we shall try to 
assess whether this approach might work in practice. 

One aspect counteracting the application of the cur- 
rent approach to the analysis of anomalous diffusion is 
that to date semiconductor devices tend to be 'too clean': 
In state-of-the-art high-mobility samples, the mean free 
path is of the same order as the low-temperature phase 
coherence length, roughly about 10 iim. In such devices, 
the phase coherent electron transport is ballistic and not 
even conventionally diffusive. We thus need to consider 
low-mobility devices, where the disorder concentration is 
increased either by doping or by lowering the separation 
between the 2DEG and the donor impurities. We expect 
that by artificially increasing the disorder, an order of 
magnitude separation between I and might be attain- 
able 1^^. Second, to observe significant deviations from 
standard diffusion, we need to be in a regime of a low 
global conductance g. In low- mobility systems showing 
integer quantum Hall transitions (when placed in strong 
perpendicular magnetic fields), the typical Coulomb en- 
ergy is low as compared to the disorder energy scale. 
In such systems, the conductance g is of order unity in 
the transition regimes, and anomalous diffusion might 
be observable by our method. Notice, however, that for 
g — 0(1), the tunneling current (like any other meso- 
scopic transport observable) will be subject to significant 
renormalization by interaction processes. More specifi- 
cally, when lowering one will run into a regime, where 
the zero-bias anomaly renormalization of the tunneling 
DoS ceases to be small. However, a quantitative analysis 
of the interplay of anomalous diffusion and interaction is 
beyond the scope of the present paper. 



IV. SPECTRAL AND PARAMETRIC 
CORRELATIONS 

In the previous section, the focus was on exploring the 
spatial structure of electronic correlation functions in ex- 
tended systems. We now turn to a complementary ap- 
plication, viz., the analysis of spectral structures in finite 
chaotic environments. Specifically, we shall consider a 
situation where one of the layers - either by top gates 
or through etching - has been converted into a quantum 
dot (QD) of characteristic size L, see Fig. ^. The 2DEG 
underneath is extended as before. 



2DEG 




FIG. 3. QD - 2DEG. 

Our objective is to analyze spectral and parametric 
correlations through the tunneling current statistics. As 
before, the bias voltage V and the parallel magnetic 
field B|| will be used to identify different contributions 
to the current-current correlator, and to detect their de- 
pendence. In addition to these control parameters, we 
will employ a perpendicular magnetic field B± to probe 
parametric correlations. 

In the finite size setup under consideration. Coulomb 
charging effects are likely to play some role. For didac- 
tical reasons, we will begin by discu ssing the idealized 
noninteracting situation in Sec. IV A. Coulomb correc- 
tions will then be considered in Sec. iVB where it will 



be argued that the impact of interactions on the current 
fluctuations sensitively depends on the parameter regime 
under consideration. 



A. Noninteracting case 

We are interested in the behavior of a chaotic quan- 
tum dot on time scales, where the electron dynamics is 
ergodic - a 'zero-dimensional' system in the standard ter- 
minology of mesoscopic physics. For a diffusive system, 
the time it takes to establish ergodic dynamics is set by 
-Erj^jJ. For nearly clean quantum dots, the ergodicity time 
depends on the specifics of the boundary scattering po- 
tential. Ergodic mesoscopic systems are tailor made to 
modeling in terms of random matrix theory (RMT) . For 
simplicity, we shall assume that the perpendicular mag- 
netic field is strong enough to globally break time reversal 
invariance. On the other hand, the field is assumed to 
be too weak to significantly affect the ballistic dynam- 
ics of the charge carriers: Rc ^ min(L, ?). Under these 
conditions, the quantum dot can be modeled in terms of 
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random matrices drawn from the Gaussian Unitary En- 
semble (GUE). Specifically, we will describe the upper 
system in terms of some A'^-dimensional Hermitian matrix 
Hamiltonian {H + where H is taken from a GUE 

and the antisymmetric matrix $ represents the magnetic 
field in a way to be discussed momentarily. The indices 
/i, v can, roughly, be interpreted as discretized spatial co- 
ordinates. (Neither the discretization in terms of N sites 
nor the specific interpretation of the indices will play any 
role throughout.) 

Within the RMT description, space-type matrix ele- 
ments will be represented as L^A(x, x') NA^j^^, and 
the integration over the coordinates of the upper sys- 
tem becomes a matrix trace, J (Px N~'^ On 
the other hand, the dynamics in the 2DEG spectrometer 
underneath is integrable-ballistic implying that, as be- 
fore, it has to be described in terms of the microscopic 
Hamiltonian of the two-dimensional electron gas. This 
type of hybrid modeling, involving RMT in combination 
with a microscopic Hamiltonian, does not pose any con- 
ceptual problems. In our basic formulas, (||), (0), and 
( p^ ) , we describe quantities assigned to the upper system 
through their RMT representations whereas the correla- 
tion functions of the lower system are given by the ballis- 
tic expression ( [T^ ) . Specifically, the RMT representation 
of the correlation function of the upper systems is 
given by 

where G;^^ = {e±iO- Ht i*/2)^J are the RMT Green 
functions, and we have assumed that the mean value 
of the perpendicular field has been absorbed into 
the Hamiltonian H. The other two correlation functions 
p[d;C] g^j.g given by similar expressions with exchanged 
coordinates. Notice that none of the functions F ac- 
tually depends on the coordinate arguments - the er- 
godic zero-dimensional nature of the dot - which is why 
no conflicts arise from the simultaneous appearance of 
RMT and truly space dependent correlation functions in 
Eq. (0). 

To be specific, let us assume - a condition easily met 
experimentally - that the electron dynamics on scales 
L is ballistic. The average tunneling conductance be- 
tween the dot and the extended system is again given 
by Eq. (||). The only difference is that, for a clean quan- 
tum dot with irregular boundary scattering, the scale T is 
given by r'^'') ^ vp/L [where the superscript '(b)' stands 
for ballistic] . 

Turning to the fluctuations of the conductance, we 
first note that for a GUE Hamiltonian, the cooperonic 
correlation function fI'-^I = 0, and we are left with 

^i/j'if' = -F''''^', where the second denotation emphasizes 
the coordinate independence. Substituting these expres- 
sions into Eq. ( [Til ) , temporarily setting the parallel mag- 
netic field to zero, and integrating over coordinates we 
obtain 



/r(b)v 

CG{V,AB)^2i — ] 



27rr(b) 



(18) 

The global suppression factor (F^'^^/eF)^ ~ (fcpi)"^ 
stems from the fact that an averaging over the area of 
the QD is intrinsic to the setup. The further suppression 
factor S/r^^'> ^ (fcpi)"^ multiplying the diffuson contri- 
bution results from the integration over the Green func- 
tions of the lower 2DEG. Unlike the density-density con- 
tribution, fI^i is weighted by Green functions taken at 
different coordinates. Integration over these arguments 
leads to the fcp-dependent suppression. 

According to Eq. (18), the fluctuations of the tunnel- 
ing conductance are linearly related to the sum of two 
ergodic correlation functions F^'^^ and - multiplied by a 
small constant - F^-^I. We next turn to the discussion 
of these objects. RMT parametric correlation functions 
have been calculated within the framework of the non- 
linear sigma model ||ll|,0 • For completeness we have in- 
cluded a brief review of these analyses in the Appendix. 
The result, an integral representation of the two functions 
pld;D] ^ given by 



oo 



-by\y-2\ _ -by{y+2)-\ 



(19) 



i^I^l is; b)-^J^ e-^y [{by\y - 2| + 1) e-^^l^'^l 



-(6y(y + 2) + l)e-^^(^+2)] . (20) 

In these equations, the correlation functions are charac- 
terized through two dimensionless parameters, s and b. 
The parameter 

s+ = 7r(U-U' + iO)/(5| 

measures the energetic mismatch between the two Green 
functions. (Throughout this section we will add to pa- 
rameters of the upper/lower system a superscript t / i: 
if necessary.) The other parameter, b ~ {NS-\)~^ tr ($^), 
describes the field difference. As shown in Ref. jl^, this 
parameter can unambiguously be related to the magnetic 
field threading the dot. Comparison of the RMT cr-model 
and its microscopic counterpart leads to the identifica- 
tion, 

26 = TT^Cct)^ , 

where C ~ <5f^((^)^) describes the sensitivity of lev- 
els to the applied field, and (j> = AB±L^ the magnetic 
flux threading the dot. For a disordered system, C is 
proportional to the conductance g. 

For general parameter conflgurations (s,6), the inte- 
grals in Eqs. (|l^) and (|2^) cannot be done in closed form. 
In a number of relevant limiting cases, however, simplifi- 
cations are possible. We first note that for zero magnetic 
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field diff'erence, A_B^ = 0, the density-density correla- 
tor F^'^\s,b = 0) is related to the two-level correlation 
function (see, e.g., Ref. |42| ) 

through 5Ri^[''l = ttS{s) - R^{s) (the index 'u' standing 
for unitary). Thus, the density-density contribution to 
the current correlation function Cj is simply proportional 
to the level number variance Y,2{V/5^). Similarly, the 
real part of the zero-field diffuson correlation function 
is given by 3?F[^1(s,6 = 0) = 7r(5(s) (where the (5-type 
dependence on s follows from the condition of particle 
number conservation) Notice that the s = singu- 
larity displayed by the two correlation functions is of no 
relevance for our theory. The reason is that for voltages 
Ay ~ (5|S < Tq- the tunneling approximation employed 
in calculating the current is no longer applicable. What 
happens for smaller voltage differences cannot be figured 
out within the framework of the present formalism. In 
principle, the (7-model description can consistently be ex- 
tended to include the effect of coherent multiple tunnel- 
ing. (A treatment neglecting the coherence of higher- 
order tunneling processes would simply lead to a complex 
energy shift, s — > s -I- 17, where 7 = ttTt /S^ and, very 
much like dephasing to a smearing of the singular low- 
energy profile of the correlation functions.) The resulting 
theory, however, would be significantly more complicated 
than the present formalism and we will not discuss it any 
further (thereby paying the price that the very low volt- 
age regime remains out of reach). 

We next focus on the dependence of the correlation 
functions on a magnetic field. There are two different 
cases to be distinguished: For small field strength, 6 <C 1, 
the effect of the magnetic field is essentially nonpertur- 
bative (in the sense that the correlation functions are 
nonanalytic in the inverse field strength ~ &~^). In this 
regime, the conductance correlation function has to be 
computed through numerical integration of Eqs. ( p^ and 
(pO|). The leading order term for small b is linear in b (i.e., 
quadratic in the field difference Ai?^). For larger fields, 
6 3> 1, a perturbative expansion in powers of b^^ and 
leads to the results 

(26-is+)2' 2b~is+ 

which would otherwise obtain from a perturbative dia- 
grammatic analysis (similar to the one outlined in the 
previous section). 

We finally ask how the two contributions F^'^l and F^^l 
contributing to the conductance correlation functions can 
be distinguished. In general, the factor (fcpi)"^ <C 1 
leads to a massive suppression of the diffuson contri- 
bution as compared to the density-density contribution. 
Since corrections of 0{{kFL)~^) have been neglected in 
applying random matrix theory, anyway, an additional 



parameter is necessary to resolve the diffuson contribu- 
tion. As discussed previously within the context of two 
extended systems, this is exactly what an in-plane mag- 
netic field does: Coupling only to the D-contribution, 
an in-plane field difference 'B^ ~ B|| — B'| can be used 

to selectively identify the F^-^l correlation function. In 
fact, for finite Bjj", the tunneling matrix elements pick up 
a phase factor that modifies the subsequent integration 
over the Green functions of the 2DEG. For small fields, 
this leads to 

Cg {V, AB^ ■ Aq) -Cg{V, AB^ ; 0) ^ 
- - (^(|Aq|L)2 ^^[^1(., b) + ((|Aq|L)4) , 

where Cg is a constant of order unity. 

The relevant field scale is one flux quantum through 
the area spanned by the linear size of the QD, L, and 
the distance between the layers, d, i.e. ^ l/{dL). 
This field scale is typically much larger than the charac- 
teristic scale of the perpendicular magnetic field B^ ~ 

Before leaving this section, let us briefiy comment on 
the connection between the present analysis and previ- 
ous experimental work by Sivan et al. | [To|] . As men- 
tioned in the introduction, Ref. |lo| investigated a setup 
dot/dot, where the second dot was very small, with ex- 
treme level quantization. (That is, in Ref. ||l^, the role 
of our 2DEG is assumed by a single quantized level of 
an ultrasmall device.) Arguing semiquantitatively, Sivan 
et al. related the statistics of the tunneling conductance 
to density-density-type parametric correlation functions. 
Indeed, the experimental data turned out to be in good 
accord with the RMT prediction, Eq. (|l9|), discussed 
above. There are two differences to the presently dis- 
cussed setup: first, the fact that in our system a 2DEG is 
used as a spectrometer implies that the second, diffuson- 
type correlation function plays a more important role 
than in Ref. ||l^]. (In fact, this type of correlation fimc- 
tion should contribute to the data of Ref. |l^, too. Due 
to the fact, that the current flows into a single level, 
however, this contribution is minute and can safely be 
neglected.) Second, one may hope that the spatial av- 
eraging involved in our formalism leads to a far-reaching 
elimination of all non-spectral structures (as opposed to 
single-level spectroscopy, where non-universal wave func- 
tion characteristics may affect the result). The price one 
has to pay is the suppression factor (fcpi)"^ that is not 
present in the dot/dot setup. 

Summarizing, we have found that the statistics of the 
dot/2DEG-conductance (in a regime of broken time re- 
versal invariance) can be described in terms of two basic 
correlation functions fI'^I and F^-^^. As compared to the 
previously discussed case of two extended systems, the in- 
formation contained in F^'^''-^^ is now purely spectral. (All 
spatial structures have equilibrated due to the ergodic- 
ity of the system.) Indeed, these two functions are fully 
universal in the sense that they depend only on the two 
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basic parameters s and b measuring bias voltage and per- 
pendicular magnetic field strength, respectively. There is 
one non-universal element affecting the conductance fluc- 
tuations, viz., a geometry-dependent factor suppressing 
the contribution of F'^I . Still the contributions of the d- 
and Z3-correlation function can be disentangled, namely, 
by measuring the dependence of the conductance on a 
parallel magnetic field. These results have been obtained 
within a noninteracting theory. We next ask what might 
change once Coulomb interactions are switched on. 



B. Coulomb interaction effects 

As far as quantum dots are concerned, the most im- 
portant manifestation of interaction phenomena is the 
Coulomb blockade: due to the repulsive interaction, it 
costs extra energy Ec to add an electron to the dot. 
For an isolated dot, this charging energy is determined 
through the self-capacitance C through Ec = e^/(2C). 
However, for a dot in close vicinity to an extended con- 
ducting system, a 2DEG, say, an excess charge on the 
dot will be compensated for by the accumulation of pos- 
itive background charge in the large system. Under such 
circumstances it is the geometric capacitance between the 
systems that determines the charging energy. This is the 
situation given presently. 

For a planar dot/2DEG setup, the geometric capaci- 
tance estimates to C = ccqL^ /d where eo ~ 8.9 • 10^^^ 
F/m and e is the dielectric constant of the filling medium. 
For GaAs, e « 10. The capacitance determines the charg- 
ing energy Ec that has to be compared with the other 
characteristic energy scales of the problem. The smallest 
scale we might hope to resolve is the single-particle-level 
spacing of the dot. (It is this energy scale on which 
non-perturbative structures in the parametric correlation 
functions are observed.) Notice that both, 6^ = (i^i^)^^ 
and Ec, scale inversely with the dot area. Thus, it must 
be the spacing between dot and 2DEG that determines 
the crossover criterion. Specifically, with i/ Ri 3 ■ 10^° 
meV~^cm~^, we find that Ec « (5| for d ~ 4 nm. Realis- 
tically, d is somewhat larger than that, i.e., of the order 
of up to a few tens of nanometers. We thus conclude 
that interaction effects are of relevance once one gets in- 
terested in low-energy structures of the order of the level 
spacing. 

In previous sections, we have described fluctuations 
of tunneling transport coefficients through disorder cor- 
relation functions J-Q{ei — £2) = (^o(ei)^o(e2)), where 
Ao(e) are the energy-dependent spectral functions and 
the subscript '0' means 'noninteracting'. (In fact, we 
have largely focused on the correlator of Green func- 
tions Fq ~ (GqGq). Presently, however, it will be 
more convenient to concentrate on the spectral functions 
themselves. The two quantities Fq and J^o are related 
through To = {AoAo)c = 25ft(Go GJ)c = 2?fiFo.) As 
mentioned above, interactions will mainly manifest them- 



selves through global charging mechanisms. This sector 
of the Coulomb interaction does not couple to the coordi- 
nate dependence of the correlation functions. To simplify 
the notation, we therefore temporarily suppress spatial 
coordinates in the notation. 

Our main goal will be to show that interaction and dis- 
order are largely separable in the analysis of correlation 
functions. Yet, unlike in the noninteracting case, it will 
no longer be sufficient to compute the zero-temperature 
correlation functions and to account for finite tempera- 
ture effects in the end, through an integration over the 
Fermi distribution functions. Instead, we have to work 
in a finite-temperature formalism from the outset. 

To model the interaction we employ the 'orthodox 
model', that is we add a charging term 

Hc = ^iNQo-Nf, (21) 

to the Hamiltonian of the dot. Here N is the preferred 
occupation number that can be set by the gate voltage. 

To incorporate this term into the model, it is con- 
venient to switch to a functional integral formulation. 
Within this approach, it is straightforward to see that 
the theory essentially splits into two sectors: one that 
describes noninteracting Green functions (albeit subject 
to some imaginary time dependent voltage u), and an 
interaction- and temperature-dependent weight function 
that controls the fluctuations of a. This approach of 
describing charging was introduced by Kamenev and 
Gefen p4| . In the following, we shall briefly review its 
main elements, and apply it to our present problem. 

Within a fermion field integral approach, the imagi- 
nary time action describing the quantum dot is given by 

J dri^idr + Ho- + Ec j dT{4,ijj - N)\ 

where tp is a time- and position-dependent Grassmann 
field and the rest of the notation is self explanatory. (Un- 
less stated otherwise, the notation ipip = J d'^ri(;{r)ip{v) 
contains a spatial integration.) Decoupling, the interac- 
tion by a Hubbard- Stratonovich transformation one ob- 
tains the effective action 

S[)p,a] = S[a] + Sd[ijj,(y], 

S[a]^ J dT{^^a^~iNa), 

Sd[i^,cr]= J drip (t) {dr + H - + ia) iP{t), 

where (t(t) is a scalar bosonic field. Next, all but 
the static component cto = P J dr a are removed from 
the fermionic action, 5d, through the gauge transfor- 
mation 4'{t) exp[— dT'{a{T') — ao)]^j{T). [That 
the static component cannot be removed has to do with 
the fact that gauging out ctq would, in general, lead to 
a violation of the time-antiperiodic boundary condition 
-0(r) = —iP{t + /?).] This makes the action of the sys- 
tem insensitive to the time-dependent components of the 
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Coulomb field. However, the (imaginary time) Green 
functions G{t), we actually wish to compute, being non- 
gauge invariant objects, pick up a gauge factor Q{t, a) 
Q{T^ao)B[T) to be specified momentarily. For tempera- 
tures larger than the level spacing, the integration over 
the static component ctq can be done in a saddle point 
approximation. As a result, ^/(t, (Tq) Qq(j,Ii) = Qoij), 
where pL has the significance of an effective (real) chemi- 
cal potential determined through the optimal occupation 
number N . 

We are thus led to consider the combination G{t) = 
Qo{t)B{t). Roughly speaking, the physics of interac- 
tions resides in the factor B{t). Disorder, the external 
fields, etc., are contained in Qo = Go{^, x'; r). This is the 
'splitting' of the zero-mode interaction theory into two 
components mentioned above. 

Transformation of Q back to frequency space ob- 
tains El 



G{ie„ 



1 



duj' de' 
'2tt2tt 



B{u')Ao{e') X 



coth ^ + tanh ^ 
icn — uj' — t' 



(22) 



where is the spectral function associated with the 
(real-time) Green function Gq and 



B[lo) 



1 



TE, 



exp[ 



ATE, 



iEl+u^)]s\nh 



2T 



This representation holds for any particular realization 
of the disorder and the external fields. 

We next turn to the discussion of the correlator of two 
spectral functions T{uo\ e) = {A{e+^)A{t - |)). Making 
use of Eq. ( p^ ) and noticing that the disorder average 
couples only to the functions Aq, we obtain 



TE, 



cosh ■ 



dW 



77) e ST 
- -I- cosh 



2T 



cosh 



~w 

T 



cosh 



(23) 



Notice that, unlike the non- interacting To depends 
only on the energy difference oj. 

Equation ( p3| ) represents the most general form of our 
result for the correlation function in the presence of a 
charging interaction. To understand the structure of this 
expression, and to identify physically distinct regimes, 
we realize that J- depends on four characteristic energy 
scales: the temperature T, the charging energy Ec, the 
bias voltage V , and an intrinsic scale eo over which the 
non-interacting correlator Tq varies. Notice that the de- 
pendence on V is implicit, through the limits imposed on 
the variables e and uj. 

To exemplify the dependence of the correlator on these 
scales, let us assume that is proportional to a 

delta function smeared over the intrinsic level broaden- 
ing Fr, i.e., To{lo) ^ (5rr(^) (as is the case, e.g., for the 



diffuson-contribution to the ergodic correlation functions 
discussed above). In this limit. 



1 



2Vm 



7(e) 



(24) 



where / is some function that decays on a scale set by T. 
This expression displays a feature common to all corre- 
lation functions T affected by charging: Formerly sharp 
energy dependences are broadened to Gaussians of width 
cr = 2\/TEc- In other words, energetically sharp features 
of J-Q are washed out and a lower limit to what can be 
resolved in an experiment is set. 

To make further progress, we have to distinguish be- 
tween two different regimes: (a) weak interaction or high 
temperature, Ec ^ T and (b) strong interaction or low 
temperature, Ec ^ T. We begin by discussing the first 
case, (a). 

For Ec <^ T, interaction corrections are small and an 
expansion to first order in the parameter Ec/T obtains 

T(uj-e)^To{Lo){l + 0{EjT)) , 

where To stands for an energy average of the noninter- 
acting correlator JFqi over a scale 2VT^ (as in Eq. (H)). 
This means that weak interactions leave the result essen- 
tially unaltered, albeit lowering its spectral resolution. 

We next turn to the discussion of case (b), Ec ^ T. 
Given that the principal setup of the theory favors low 
temperatures, this regime is more relevant than (a). On 
the other hand, it also has to be kept in mind that the 
applicability of the theory is limited to temperatures 
T > 6. Thus, the structures discussed below apply to a 
temperature regime S <T <^ Ec- 

In the extreme limit T ^ [Q, the interactions pro- 
duce a hard Coulomb gap, i.e. T = To 9{e — \uj\/2 — Ec). 
At finite temperatures the gap becomes softer but its 
essential features remain robust: For small bias V, the 
correlation function is strongly suppressed while at large 
bias there are only s mall changes. The relevant limits are 
(bl) V < Ty^EjT and (b2) T<s:Ec<^V. 

In the first case, when the applied voltage is small, 
Eq. ( p3| ) simplifies to 

Tiuj;e) ~ cJPo(0)e-^ (cosh |; + cosh ^), 

i.e., the correlator is exponentially suppressed. Here c is 
a factor depending algebraically on T/ Ec. 

In the opposite case, when the applied voltage is large, 
one obtains 

T{u);e) ~ To{uj) - 

-e^^~T^ [ToiuJ-] e~5T + To{'jJ+) e^) , 

where the dominant contribution To{iij) is the noninter- 
acting To smeared out over energies 2^/TEc as before 
and ijj± — UJ ± 2Er. 
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Summarizing, large charging energies, > T, change 
the noninteracting theory in two different ways. First, to 
avoid the Coulomb blockade, large bias voltages V > Ec 
have to be applied. Second, even for those voltages, a 
lower limit 2^/TEc on the maximal energetic resolution 
of the theory is imposed. For the reasons outlined above, 
we believe that, at least for sufficiently weak magnetic 
fields and impure samples, the effect of other interaction 
mechanisms is relatively minor. At any rate, since the 
functional dependence of the Coulomb blockade correc- 
tions follows from Eq. ( p3| ) , the applicability of the theory 
can be put to test. 



V. CONCLUSIONS 



In this paper, we have introduced an approach to ex- 
ploring the three basic two-particle correlation functions 
of mesoscopic physics, the generalized diffuson F^-^\ the 
cooperon fI'-^I, and the density-density correlator F^'^\ 
The basic idea is to monitor the current flowing between 
two parallel 2DEGs as a function of three qualitatively 
different control parameters - a parallel magnetic field, a 
perpendicular field, and a bias voltage. 

As compared to standard transport measurement ar- 
chitectures, the most important advantage of the ap- 
proach is that electronic correlations are detected with- 
out disturbing the 'bulk' system through local contacts. 
Instead the entire planar electron system acts as an 
'extended contact' whose spatial structure is, nondis- 
turbingly, scanned by means of the two magnetic fields; 
spectral electronic structures are resolved by measuring 
the bias voltage dependence of the current. Importantly, 
the relevant information carried by the tunneling current 
is solely contained in its fluctuations. For example, we 
have shown that the Fourier transform of the current- 
current correlation function with respect to the external 
field, directly obtains the two spatially resolved trans- 
port functions F^^''^^ In contrast, previous analyses of 
magnetotunneling currents, both experimentally and the- 
oretically, focused on the average current profile that is 
unrelated to any 'mesoscopic' type of information. 

To exemplify the usefulness of the approach, we have 
considered two different applications: tunneling between 
two extended 2DEGs and tunneling from a quantum dot 
geometry into a 2DEG. As for the former, we have shown 
how, at least in principle, the exponents characterizing 
localization/delocalization transitions can be extracted 
from the current data. In contrast, for a quantum dot 
with ergodic dynamics, the focus is on spectral rather 
than on spatial structures. We have explicitly worked 
out the connection between the tunneling current statis- 
tics and various parametric correlation functions (a con- 
nection previously used on a semiquantitative basis to in- 
terpret the data of the experimentreported in Ref. |10|) 
and discussed how these correlations can be monitored 
by changing external field and bias voltage. 



As with any other tunneling setup. Coulomb interac- 
tions are likely to change the outcome of the noninter- 
acting theory. We have provided evidence in favor of the 
Coulomb blockade being the most relevant interaction 
mechanism. The presence of the Coulomb blockade will 
result in two principal effects: first, it forces one to use 
bias voltages in excess of the blockade threshold. Sec- 
ond, the spectral structure of the correlation functions is 
washed out. This reduces the information content that 
can be extracted from the tunneling current statistics. 
The extent to which these obstructive mechanisms affect 
the theory is set by the Coulomb charging energy. In the 
present system, the latter is largely determined by the 
inverse of the interlayer capacitance which, owing to the 
extended geometry of the systems, is small. Thus, charg- 
ing phenomena will be less pronounced than in small is- 
lands. All in all, we believe that for realistic values of 
the relevant system parameters, a significant parameter 
range over which electronic correlations can be resolved 
through the current approach remains. 
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APPENDIX A: RMT CORRELATION 
FUNCTIONS 



In this appendix, we briefiy recapitulate the essential 
steps of the the nonlinear a model (NLcrM) analysis of 
RMT parametric correlation functions. The goal is to 
compute correlation functions of an advanced and a re- 
tarded Green function taken at different energies and 
magnetic fields. To model the field difference, we assume 
that the Hamiltonian of the retarded/advanced Green 
function is given by H ±i^/2. Here $ is some antisym- 
metric matrix that can be modeled in a variety of dif- 
ferent ways We are thus led to consider correlation 
functions of the type 



2 



H-i^/2) 
^ 2 



1 



i^/2) 



Correlation functions of this type can conveniently be 
calculated from the functional integral 



(Al) 



S = i^K'/^ie - iJ - - i$)A}Ai/V . 



Here A = a^^ is a Pauli matrix in the 'advanced-retarded' 
space, u)^ = uj + iQ, and the ellipses stand for pre- 
exponential terms specific to the index configuration of 



13 



the correlation function under consideration Et],^ ■ Fol- 
lowing a by now standard procedure |^ , ^ , pl| 7the aver- 
age over the RMT Hamiltonian then leads to the effective 
action of the NLcrM, 



S-eff = ^ STr(QA) - ^ STr(QA)2 , 



(A2) 



where Q is a four-dimensional supermatrix subject to 
the constraint = 1. The parameters, s and b, appear- 
ing in this expression measure the 'mismatch' of our two 
Green functions. Their relation to system parameters is 



discussed in Sec. [TV A 

At this stage the computation of the parametric 
correlation functions fI'^'-^''^! amounts to an integral 
over the four-dimcnsional matrix Q. The explicit 
calculation leads to F^'^^ = and the two 

one-dimensional integral representations displayed in 
Eqs. (|l|) and (f^) in the main text. 
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